Phase Transition in Potts Model with Invisible States 
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We study phase transition in the ferromagnetic Potts model with invisible states that are 
(3JQ! added as redundant states by mean-field calculation and Monte Carlo simulation. Invisible 

states affect the entropy and free energy, although they do not contribute to the internal 
^* ■ energy. A second-order phase transition takes place at finite temperature in the standard 

g-state ferromagnetic Potts model on two-dimensional lattice for q — 2, 3, and 4. However, 
our present model on two-dimensional lattice undergoes a first-order phase transition with 
spontaneous (/-fold symmetry breaking (q — 2, 3, and 4) due to entropy effect of invisible 
states. The model is fundamental for the analysis of a first-order phase transition with 
spontaneous discrete symmetry breaking. 






It has been a central issue in statistical physics to study the nature of phase 

transition with spontaneous breaking (SB) of symmetry since the establishment of 

the Ginzburg-Landau theory. Recently, the concept of phase transition has not 

just been for natural phenomena but also for problems on information science, for 

c/2 , example, coloring problem and the so-called satisfiability problem. Study on phase 

transition with SB of symmetry has become increasingly important. The Ising model 

rH . and Potts model have been established as minimal models for the analysis of phase 

*^ \ transition with SB of discrete symmetry. The spin- 1/2 ferromagnetic Ising model 

'"O ' without an external magnetic field undergoes a second-order phase transition with 

SB of twofold symmetry at finite temperature. Similarly, in the ferromagnetic q- 

q \ state Potts model with no external magnetic field, g-fold symmetry is spontaneously 

broken at finite temperature. It is important to investigate the relation between the 

£SJ . order of the phase transition and the symmetry that breaks spontaneously at the 

transition point. 
_j ■ In the g-state Potts model, all spins take q states and the interaction between 

spins is represented by Kronecker's delta. 1 ' Since the Potts model for q = 2 is equiv- 
alent to the Ising spin systems, it is regarded as the straightforward extension of the 
I/"") ' Ising model. In the ferromagnetic g-state Potts model on two-dimensional lattice, a 

second-order phase transition occurs at finite temperature when q < 4, while a first- 
order phase transition occurs when q > 4.^ The Potts model has been studied by an- 
alytical approaches as well as numerical methodsPSiSaJaJSJIDnSilSJSIJSJOBJJaEDJaDJTaiSIJTaiaj 

To analyze the phase transition with SB of g-fold symmetry, it is enough to consider 
rS ' the g-state ferromagnetic Potts model in many circumstances. The Potts model has 

C3 ■ succeeded in the analysis of phase transitions in experimental systems as well as 

theoretical systems™ Recently, a number of counterexamples have been found in 



in 



o 



*' E-mail:r. tamura@issp.u-tokyo.ac.jp 

**' E-mail:shu-t@alice. math. kindai.ac.jp 

***) E-mail:kawashima@issp.u-tokyo. ac.jp 



2 Letters 

some kinds of two-dimensional modelP^^'^* Even if systems undergo a phase 
transition with SB of (/-fold symmetry, the order of the phase transition of these sys- 
tems is different from that of the ferromagnetic ^-state Potts model. To study these 
phenomena in a more systematic way, we introduce redundant states in the standard 
(/-state Potts model. Now, we assume that redundant states do not contribute to 
the internal energy, while they affect the entropy and free energy. In this paper, we 
focus on the consequence of introducing redundant states on the nature of the phase 
transition. 

There are two main types of model where the order of the phase transition 
changes by controlling parameters. Typical examples of the first type are the Blume- 
Capel modeP^'I^Sl 1 and Blume-Emery-Grimths (BEG) modell^ In these models, 
the order of the phase transition changes when the parameter set such as a chemical 
potential is manipulated. The other type is the model where the order of the phase 
transition varies if the number of microscopic states changes (e.g., the Wajnflasz 
model) j22}J32Iil2D The Wajnflasz model is regarded as the generalized Ising model 
in which the numbers of +l-states and — 1-states are different. In this model, the 
entropy-driven magnetic field is induced due to the bias of the numbers of + 1-states 
and —1-states. As a result, the order of phase transition can be changed. However, 
this phase transition does not accompany the breaking of Z<i symmetry. Although 
our present model is similar to the latter case, the phase transition accompanies the 
breaking of (/-fold symmetry. Thus, our model is specialized for the analysis of the 
relation between the SB of symmetry and the order of phase transition. Since our 
present model is the generalized Potts model as we discuss later, it can be applied 
to the analysis of SB of (/-fold symmetry. 

Our model is defined by the following Hamiltonian, 



H = ~ J ^2 Ss i> s i 5Z K,aS aj ,a, (1) 

(i,j) a=l 

Si = !,■■■ ,q,q + l,-~ ,Q + r, (2) 

where 5 Si , s . denotes Kronecker's delta. We consider the case of ferromagnetic inter- 
action (J > 0). Hereafter, we call this model ((/,r)-state Potts model, where N is the 
number of spins. The total number of microscopic states of this model is (q + r) , 
whereas that of the standard (/-state Potts model is q N . We call redundant states, 
q+1 < Si < q + r, "invisible states". Following Eq. (pfj) , if and only if 1 < Sj = Sj < q, 
the interaction is nonzero. This Hamiltonian for r = is equivalent to the standard 
(/-state Potts model. Note that the internal energy does not change when the invis- 
ible states are introduced. Thus, the ground state of this system has no invisible 
state and the number of ground state is q. It is clear that a phase transition with 
SB of (/-fold symmetry occurs in the present model. It is easy to show that the (q,r)- 
state Potts model is equivalent to the (g,l)-state Potts model with an external field 
depending on the temperature. Introducing new spin variables, Oi = 0, 1,2, • • • ,q 
related to the original spins by Oi = s, (if Sj = 1, •••,(/), o% = (otherwise) as the 
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invisible state, and taking the trace over Sj, we obtain the following Hamiltonian: 

q 
W = -j'YlSffuffj J~]8er i ,aScr i ,a ~ Tlog r ^ 5 azfi , (3) 

(i,j) a=l i 

where T represents temperature and we set the Boltzmann constant k-g to unity. 
Note that the partition function of Eq. (pQ) and that of Eq. ([3]) are the same. The 
total number of microscopic states of this representation is (q + 1)^. The second 
term in Eq. ([3]) comes from reducing the degrees of freedom of invisible states in 
the original Hamiltonian given by Eq. (pQ), and it corresponds to the entropy effect. 
This term represents the chemical potential for the invisible states, which induces 
the creation of invisible states at high temperature. 

In this paper, we consider the order of the phase transition of the ferromagnetic 
(g,r)-state Potts model on a square lattice with periodic boundary condition. We 
focus on the cases of q = 2, 3, and 4, because the standard ferromagnetic Potts model 
on two-dimensional lattice has a second-order phase transition for these values. 

Before considering the nature of the phase transition on a two-dimensional lat- 
tice, we analyze the phase transition of this model by mean-field calculation based 
on Bragg- Williams approximation with two parameters that represent the order pa- 
rameter and density of invisible states. We first consider the (2,r)-state Potts model 
that can be mapped on the BEG modelP3 

^BEG = -~ J2 {tit, + tjt]) -Dj2{l~tj), (4) 

i, = +1,0,-1, (5) 

where D represents the crystal field. Hereafter, we take J as the energy unit. When 
D = Tlogr, this Hamiltonian is equivalent to the (2,r)-state Potts model. Now, we 
consider the case that the coordination number is four corresponding to a square 
lattice. The left panel of Fig. [1] shows the phase diagram of the BEG model obtained 
by mean-field calculation. The dotted curve and solid curve indicate first-order and 
second-order phase transition points, respectively. Thus, the decrease in temperature 
in the (2,r)-state Potts model corresponds to a tilting trajectory followed by D = 
T log r in the phase diagram shown in the left panel of Fig. [TJ As r increases, the 
transition temperature decreases. In the right panel of Fig. [H we show the latent 
heat and transition temperature as functions of r for q = 2, 3, and 4. 

A second-order phase transition takes place only if (q,r) = (2,1), (2,2), and (2,3). 
Note that the standard g-state ferromagnetic Potts model (q > 3 and r = 0) has a 
first-order phase transition as far as mean-field calculation.^ As r increases, the 
latent heat increases and the transition temperature decreases. For r — > oo, the 
latent heat and transition temperature approach 2 and 0, respectively. From the 
mean-field calculation, we expect that the order of the phase transition changes by 
introducing invisible states. 

We also investigate the nature of the phase transition of the (<7,r)-state Potts 
model by Monte Carlo simulation using the standard Metropolis method. Each run 
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r=2 — r=4 



r=io; 



r=3 



r=10 -- r=10 J 




Fig. 1. (color online) Results of the mean- field calculation for the BEG model (Eq. Q). (left) 
Phase diagram. The dotted curve and solid curve represent first-order and second-order phase 
transition points, respectively. Tilted lines represent the trajectories of decreasing temperature 
for r — 2,3,4, 10, 10 2 , and 10 3 (from bottom to top), (right) Latent heat AE and transition 
temperature T c as a function of r. 
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Fig. 2. (color online) Specific heat C and the square of the order parameter \m\ as functions of 
temperature of the (3,25)-state Potts model. The error bars are omitted for clarity since they 
are smaller than the symbol size. 



contains 10 7 — 10 8 Monte Carlo steps per spin at each temperature. We make 8-16 
independent runs for each size to evaluate the statistical errors. Figure [2] shows the 
specific heat C and square of the order parameter of the (3,25)-state Potts model as 
a function of temperature. 
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Fig. 3. (color online) Probability distribution of the internal energy of the (3,25)-state Potts model. 
The error bars are omitted for clarity since they are smaller than the symbol size. 



It is convenient to introduce another representation of Kronecker's delta as 
$a,p = [l + (q — 1) e a ■ e@~\ /q, where e a (a = 1, 2, • • • , q) represents q unit vec- 
tors pointing in the q symmetric direction of a hypertetrahedron in q — 1 dimen- 
sions™ The conventional order parameter for the g-state Potts model is given by 
m = Y2i=i e -° i 1^ ■ Note that in the (g,r)-state Potts model, we adopt the standard 
order parameter and the relation that e°" i=0 = for the invisible state. As the system 
size N(= L 2 ) increases, the maximum value of the specific heat increases linearly. 
The square of the order parameter \m\ , which detects SB of (/-fold symmetry, takes 
a finite value. These facts indicate that a first-order phase transition with SB of 
g-fold symmetry occurs at finite temperature. 

To obtain the transition temperature and latent heat, we calculate the proba- 
bility distribution of the internal energy and apply finite-size scaling for a first-order 
phase transition following the method in Ref . 36) . Figure [3] shows the probability 
distribution of the internal energy P (E), which has a form of bimodal distribution. 
The bimodal distribution of energy histogram is considered to be a typical behavior 
of a first-order phase transition. The maximum values of the size-dependent specific 
heat C max (L) and the transition temperature T C (L) depending on the lattice size 
defined by C(T C (L)) = C max (L) are obtained by the reweighting method. We adopt 
the finite-size scaling relations for a first-order phase transition in a d-dimensional 
systemPSJ 



T C (L) = 



aL- d + T c , 
{AE) 2 L d 
4T? ' 



(6) 
(7) 



where AE and T c represent the latent heat and transition temperature for infinite 
systems, respectively. The top panel of Fig. 2] shows T C (L) as a function of L~ 2 . We 
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Fig. 4. (color online) Finite-size scaling for the (3,25)-state Potts model. (Top panel) Transition 
temperature depending on the lattice size as a function of L~~ 2 . (Bottom panel) Maximum 
value of the specific heat as a function of L 2 . Dotted lines represent fitting curves using data of 
L > 48. The error bars are omitted for clarity since they are smaller than the symbol size. 



fit the data using Eq. ([6]) and obtain the transition temperature T c = 0.59630(1). 
This value is lower than that of the standard 3-state Potts model. 2 ' The bottom 
panel of Fig. H] shows C max (L) as a function of L 2 . We fit the data using Eq. 
and we obtain AE = 0.81(2) using the transition temperature given by the above 
calculation. Using the same technique, we obtain the transition temperature and 
latent heat of (2,30)- and (4,20)-state Potts models which are shown in Tabled We 
confirmed a first-order phase transition for these cases as well. 

Table 1. Transition temperature and latent heat for (2,30)-, (3,25)-, and (4,20)-state Potts models. 



(q,r) 



AE 



(2,30) (3,25) (4,20) 



0.57837(1) 0.59630(1) 0.61683(1) 
1.02(2) 0.81(2) 0.68(2) 



The (q,r)-state Potts model is expected to be a fundamental model for analyzing 
a first-order phase transition with SB of q-fold (q = 2,3, and 4) symmetry in a two- 
dimensional lattice. It is found that a first-order phase transition occurs at finite 
temperature on two-dimensional fully frustrated continuous spin systems with SB 
of threefold symmetry I2aJ 42ZJ 43S> Since these models are based on continuous spin 
systems, these have additional degrees of freedom in addition to three states. Thus, 
it is possible that these degrees of freedom can be regarded as the existence of intrinsic 
invisible states. Furthermore, since the Potts model is applied to a broad range of 
science, e.g., information science l2ZJil2§IJ32I we a i so expect that the (g,r)-state Potts 
model is widely applicable for information science and technology. 

To summarize, we have studied a finite-temperature phase transition in the 
ferromagnetic Potts model with invisible states on two-dimensional square lattice. 
The invisible states are introduced as redundant states that affect the entropy and 
free energy. However, they do not contribute to the internal energy. We found 
that a thermal-driven first-order phase transition occurs in this system for q = 2, 
3, and 4 due to the entropy effect of invisible states. We also found that the latent 
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heat increases with the number of invisible states r in the mean-field calculation 
as stated above. Thus, we expect that a first-order phase transition occurs for 
larger r than those used in the Monte Carlo simulation. To find a boundary of r 
between occurring a second-order and a first-order phase transition is an issue in 
the future. The present model can be regarded as the generalized Potts model that 
is fundamental for analysis of spontaneous discrete symmetry breaking. We believe 
that the present model provides novel insight into such kind of phase transition. 
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